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ABSTRACT 

The results of an optical spectroscopic survey of a sample of young stellar objects 
(YSOs) and pre-main sequence (PMS) stars in the Lupus Clouds are presented. 92 
objects were observed with VLT/FLAMES. All of those objects show IR excess as 
discovered by the Spitzer Legacy Program "From Molecular Cores to Planet-Forming 
Disks" (c2d). After reduction, 54 spectra with good signal-to-noise ratio are spectrally 
classified. Effective temperatures and luminosities are derived for these objects, and 
used to construct H-R diagrams for the population. The sample consists mostly of 
M-type stars, with 10% K-type stars. Individual ages and masses are inferred for 
the objects according to theoretical evolutionary models. The mean population age is 
found to be between 3.6 and 4.4 Myr, depending on the model, while the mean mass is 
found to be ^0.3 Mq for either model. Together with literature data, the distribution 
of spectral types is found to be similar to that in Chamaeleon I and IC348. The Ha 
line in emission, found in 49% of the sample, is used to distinguish between classical 
and weak-line T Tauri stars. 56% of the objects show Ha in emission and are accreting 
T Tauri stars. Mass accretion rates between 10~® and 10~^^ Moyr~^ are determined 
from the full width at 10% of the Ha peak intensity. These mass accretion rates are, 
within a large scatter, consistent with the Mac relation found in the literature. 
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1 INTRODUCTION 

Lupus Molecular Clouds is the generic denomination of a 
loosely connected concentration of dark clouds and low-mass 
pre-main sequence stars located in the Scorpius-Centaurus 
OB association at 16''20'" < a < 15''30™ and -43° < 
5 < —33°. Due to its large size, close distance (d = 150- 
200 pc) and substantial mass of molecular gas, the Lu- 
pus Clouds have been subj ect of many studies at all wave- 
lengt hs o ver t he years (e.g. iHughes et al] 1994 : Men'n et al.l 



20091 : 



20081: see ICom cron 2008j for a review; IComeron et al. 
Tothill et alTBooal . 

Lupus is one of the five clouds selected by the Spitzer 
Legacy Program "From Molecular Cores to Planet- Forming 
Disks", also referred to as c2d (Cores to Disks, Evans et 
al. 2003). Using 400 hours of observations and all three in- 
struments of Spitzer, the c2d program studies the process of 
star and planet formation from the earliest stages of molec- 
ular cores to the epoch of planet-forming disks. The five 
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observed clouds cover a range of cloud types broad enough 
to study all modes of low-mass star formation, and large 
enough to allow statistical conclusions. A rich population of 
low-mass young stellar objects (YSOs) still surrounded by 
their circumstell ar material has been dis covered by Spitzer 
in these clouds llHarvev et al.l l2007al lbl: iMerm et all I2OO8I : 
lEvans et all 120091 '). 

Disks around YSOs, called protoplanetary disks, consist 
mainly of dust and gas. Dust particles of sub-^m size domi- 
nate the disk opacity, making it easily observable at infrared 
and (sub-)millimeter wavelengths by re-emitting some of the 
received stellar radiation. Dust emission is also temperature 
dependent with colder dust emitting at longer wavelengths 
than warm dust. Protoplanetary disks evolve in time, ending 
up in different scenarios like a planetary system or a debris 
disk, that may or may not also harbor planets. 

The stellar radiation spectrum depends on the effec- 
tive temperature of the star. Because of the properties of 
dust emission, the IR excess in a spectral energy distribu- 
tion (SED) of this star-|-disk system provides information 
about the geometry and properties of the dusty disk. To 
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separate star and disk emission, the stellar characteristics 
need to be known. The inner regions of the circumstellar 
disk are disrupted by the stellar magnetosphere. The stellar 
magnetic field creates channels through which material can 
flow from the disk onto the star. This produces atomic lines 
from hydrogen, calcium, oxygen, etc. The strongest emission 
line is Ha at 6562.8 A. 

The Lupus complex is one of the largest low-mass star 
forming regions on the sky, with four main star forming sites, 
referred to as Lupus LH-IH-IV. Furthermore, Lupus III con- 
tains one of th e richest associations of T Tauri stars (see 
IComer6n|[200g for a review). The work presented here con- 
centrates on data from Lupus I, III and IV because those 
were the regions observed by the c2d program. Their star 
formation rates ar e 4.3, 31 and 4.5 Mq Myr"'^, respectively 
l|Evans et al.ll2009l l. Distances of 150 ±20 pc for Lupus I and 
IV and 200 ± 2 pc for Lupus III are assumed in this work 
l|Comer6nll2008l ). 

We present an optical spectroscopic survey designed to 
characterize the young stellar population of Lupus I, III and 
IV, as observed by the c2d program. A similar study on a 
subsample is performed by Alcala ct al. (2011). Their data 
from an X-shooter survey are still preliminary, but suitable 
for a parallel study to ours. Section [2] describes the selec- 
tion criteria for the sample and Section [3] the observations 
and data reduction. Using libraries of standards from the 
literature, spectral types and effective temperatures are ob- 
tained in Section 13] Combined with additional photometric 
data, the SEDs of the objects can be built. In Section [5] the 
stellar and disk luminosities are calculated from the SEDs. 
In Section [S] the stars are placed in Hertzsprung- Russell di- 
agrams and individual masses and ages are derived based 
on theoretical tracks. The mass accretion rates, determined 
through the Ha line, are presented in Section [T] In Section 
[8] the young stellar population of Lupus is discussed in con- 
text with other regions. Finally, in Section|§l the conclusions 
from this work are stated. 



2 SAMPLE SELECTION 



iMerfn et al.l l|2008l ) used the Spitzer c2d point source catalog 
to identify the YSO population in the Lupus clouds. Objects 
are classified as YSO if they show an IR excess in the SED. 
To obtain an optimal separation between youn g stars, back- 
groun d galaxies and Galactic post-AGB stars, iMen'n et al.l 
(1200a) used the selection criteria de veloped b y t he c2d team 
on its official point source catalog (Evans et al. 20071 ). The 
method relies on an empirical probability function that de- 
pends on the relative position of any given source in sev- 
eral color-color and color-magnitude diagrams where diffuse 

boundaries have been deter mined. 

The lis t of YSOs f r om lEvans et al.1 l|2007h was then 
adapted bv lMerfn et al.l (|2008l ). Visual inspection was per- 
formed to subtract suspected galaxies or binaries, l eaving 
the list with 94 YSOs. The final fist of lMerm etlo] (|2008i ) 
was merged with 65 pre-main sequence (PMS) stars and 
PMS candidates. Here, the term PMS star is used for other 
objects added to the list whose youth had already been con- 
firmed using other observational techniques, mainly opti- 
cal spectroscopy. If an object has not been spectroscopically 
confirmed as young but it was selected by its optical and 



near-IR photometry as such, it is labeled as a PMS can- 
didate. This final list of 159 young objects is used for the 
observations presented here. 



3 OBSERVATIONS AND DATA REDUCTION 

The data presented here were taken in the second half of the 
nights of 20 - 25 February 2008 with the Very Large Tele- 
scope (VLT) and the instrument FLAMES/GIRAFFE (ID: 
080C.0473-A, PI: Oliveira). The instrument was used in the 
MEDUSA mode, with wavelength coverage of 6437 - 7183 A, 
and spectral resolution of 0.79 A. This wavelength range was 
chosen for containing temperature sensitive features, useful 
for spectral classification. Additionally, it covers the Ha line, 
an accretion diagnostic. MEDUSA has 135 fibers available, 
each with an aperture of 1.2". In total, 250 stars in 19 fields 
were observed (on average: 14 stars per field). Of those, 158 
are field stars, leaving the Lupus science sample with 92 ob- 
jects. An overview of the observations can be found in Table 
[1] The different exposure times are adjusted to the mean 
magnitude of the objects in a given field to avoid saturation 
of the brightest sources. 

The VLT data pipeline, GASGANO, was used for 
the data reduction. For each observation night, a set of 
5 dark frames, 3 fiat-fields and 1 arc frame are pro- 
duced. GASGANO performs bias subtraction, flat-fielding 
and wavelength calibration. Flux calibration was not re- 
quired for the science goals since the bands used for classi- 
fication are close to each other. Further spectral extraction 
was performed within IDL. In each observed field, unused 
fibers were placed in random "empty" sky positions. Those 
positions are not really "empty" sky, since the cloud itself 
can contribute to the sky with a lot of emission lines. In each 
field, all the sky spectra are combined within IDL into one 
master sky spectrum, with a 2(j clipping. This combined sky 
is then subtracted from the science spectra to obtain spectra 
that are ready for analysis. 

The spectra are sorted into four categories: good spec- 
tra, non-detection with a mean fiux around zero, featureless 
spectra and not useful spectra (marked as G, U, F and O 
respectively in Table [1} . The two spectra classified as not 
useful for our science purpose are objects # 14 and 53. 
Object 14 was discovered to be a galaxy (see section I4.2|l 
while object 53 was observed with a broken fibre, resulting in 
an unclassifiable spectrum. The non-detections (marked 'U', 
23 objects) and featureless (marked 'F', 13 objects) spectra 
could also not be spectrally classified. This leads to 54 good 
classifiable spectra, out of which 26 have a signal-to-noise 
ratio greater than 20. 



4 SPECTRAL CLASSIFICATION 
4.1 Method 

From the obtained spectra, the YSOs were classified by com- 
paring their spectral features with libraries of standard stars. 
To match the high resoluti on of th e FLA MES spectra, high 
resolution standards from iMonted l|l998h were used. How- 
ever, standards for all spectral types are not available in 
this library, so that only a range of spectral types can be 
derived from it. To further the analysis and narrow down 
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the spectral type, a low resolution (3.58 A) library was used 
as well (G. Herczeg, private communication, 2010). Once a 
spectral type range has been determined from the high res- 
olution standards, the low resolution library is used for a 
finer determination due to the availability of standards of 
nearly all spectral types in this library. 

The most prominent features in late type stars are the 
Titanium Oxide (TiO) absorption lines at 7050 - 7150 A. 
Mostly K-type and M-type stars have these features because 
they are cold enough for TiO to exist. The TiO feature at 
these wavelengths shows three absorption bands, which are 
deeper for colder objects. The Li I line in absorption at 6707 
A is also a common feature in young stars. This line was 
detected in some objects (see Table [2]), but was not used to 
spectrally classify the objects. 

To compare with the high resolution standards, the sci- 
ence spectra are normalized (divided by the continuum) and 
overplotted on the normalized standard spectra of different 
spectral types (see Figure [T]). A spectral range is visually de- 
termined by the best match for the TiO bands. The method 
for the low resolution standards is somewhat different be- 
cause that library is not normalized. First, the resolution 
of the science spectra is lowered to that of the standard 
(3.59A). This is done by convolving the spectrum with a 
Gaussian profile. The comparison of the spectra with the 
low resolution standards then happens in two ways. One 
way consists of scaling the spectra to the standards by an- 
choring their fluxes at certain wavelengths (6500, 7020 and 
7050 A). An example is shown in Figure[2](top left, top right 
and bottom left). In the alternative method a scale factor 
is determined by first dividing both spectra and then tak- 
ing the mean of those values, ignoring extreme features (like 
Ha emission lines) in the process. This factor is then used 
to scale the original spectra to the standards (see bottom 
right panel of Figure [2|. Both methods agree very well, pro- 
ducing four different plots for each spectral type within the 
determined range. The correct spectral type is again visually 
determined by the best match at the TiO bands. The typical 
uncertainty in the spectral classification is one sub-class. 



4.2 Special spectra 

Some spectra present special features besides the temper- 
ature sensitive ones used for spectral classification. These 
objects were inspected more closely and are differentiated 
in three special types of spectra: 

• Ten spectra show [He l] emission lines at 6677.6 and 
7064.6 A (objects # 7, 52, 76, 82, 83, 84, 85, 87, 91, 93 - see 
examples in Figure [S]). All of these objects also show strong 
Ha in emission. Strongly accreting objects produce emission 
lines other than Hq, such as this line. It is concluded that 
these lines are a sign of accretion, as is Hq (see Section [7] 
for more information on this subject). 

• Four spectra (objects #50, 59, 78 and 88) show the 
[Nil] doublet at 6548.4 and 6583.4 A and the [Sii] doublet at 
6715.8 and 6729.8 A (see examples in Figure These lines 
are indicative of ionised gas. Par-Lup3-4 (#50) and Szl02 
(#78) are related to Herbig-Haro (HH) objec ts HH600 and 
HH228 respectively (|Wang fc Henning||2009l '). They classi- 
fied these objects to be HHs with the [Sii] doublet line at 
6717/6731 A. For Par-Lup3-4 only these emission lines were 
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Figure 1. Example of the classification method of SST-Lup3-1 
with high resolution standards. The black curves represent the 
high resolution standard models and the orange curves the sci- 
ence spectrum. At the left of each standard its spectral type is 
indicated. The range of plausible spectral types is determined to 
be M4 - M8.5 
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Figure 2. Example of the classification methods of SST-Lup3- 
1 with the low resolution standard of a M6 main sequence star. 
The first three panels show the first method where the spectra 
are anchored at different wavelengths. The last panel shows the 
second method. 



detected. This is consistent with iHuelamo etU] (I2OI0D who 
show that this source has a close to edge-on disk. This would 
make the source appear underluminous. 

The two other objects that show similar spectra, 
SSTc2dJ160708.6-391407 (#59) and Szl33 (#88), are not 
yet classified as Herbig-Haro, and are thus designated as 
Herbig-Haro candidates. 

• One spectrum, IRAS15589-4132 (#14), shows lines at 
6757.95, 6774, 6794.8, 6931.5 and 6947.8 A, and no other 
features (see Figure [5]). These lines are identified as the 
doublets [N 11] (at 6548.4 and 6583.4 A) and [S 11] (at 
6715.8 and 6729.8 A) and the strongest Ho line (at 6562.8 
A) at a redshift z = 0.032. This suggests that the object 
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Figure 3. Sample of spectra (#82 on left panel and #7 on right 
panel) with strong Ha (6562.8 A), and He I (6677.6 and 7064.6 
A) emission lines. 
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Figure 4. Sample of spectra (#78 on loft panel and #59 on 
right panel) with [Nil] and [Sll] doublets (at 6548.4 and 6583.4, 
and 6715.8 and 6729.8 A respectively). Those sources may be 
associated with Herbig-Haro objects. 



is a red-shifted galaxy rather than a YSO in the Lupus 
Cloud. In the NASA/IPAC Extragalactic Database (NED), 
there is a galaxy classified at the position of this object: 
2MASSJ16022165-4140536. For further analysis, this object 
will no longer be considered part of the sample, leaving 91 
sources in total. 




Spectral Type 

Figure 6. Distribution of spectral types found in this work. The 
distribution peaks at M4-M6. 



4.3 Spectral Types 

Following the method described above, the resulting spec- 
tral types with their error ranges are listed in Table [2] The 
objects are mostly M-type stars, but 10% are K-type. The 
histogram in Figure [6] shows the spectral type distribution, 

peaking at M6. 

40 of these 54 objects were designated bv lMerm et al.l 

l|2008l ) as YSOs. Spectral classification is found in the liter- 
ature for 31 of the 54 objects (57%). The literature classifi- 
cations match well (within one sub-class) with those derived 
he re. One object, however, differs: object # 26 is classified 
bv lLopez Martf et al.] (|2005l ) as M5, but as KO (range from 
KO to K2) here. Their spectral type determination of M5 
was based on photometry, which is not as reliable as the 
spectral classification performed here. For further analysis, 
the values derived in this work are used. 
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Figure 5. Spectrum of IRAS15589-4132. The emission lines cor- 
respond with the spectral lines [N o] (at 6548.4 and 6583.4 A) 
and [S n] (at 6715.8 and 6729.8 A). This object is identified as a 
galaxy at redshift z = 0.0324. 



4.4 Effective Temperature 

Because spectral types are directly related to effective tem- 
perature T^ff, the temperatures of the stars classified here 
can be obtained (see TablefSjl . For spec tral types in the range 
K - MO, the relationship adopted bv iKenvon fc Hartmann 
( 1995) was used, while the relationship from lLuhman et al. 
(2003) was adopted for types later than MO. The uncertain- 
ties in the temperature are determined directly from the 
spectral type uncertainties. 



5 SPECTRAL ENERGY DISTRIBUTIONS 

For a given star with known spectral type (and therefore 
Tf.ff), an SED can be constructed using the correct stel- 
lar atmosphere model and additional photometric data. The 
templates for stellar atmosph eres used here are the NextGen 
Model Atmospheres (Ifauschildt et all 1 19991 : lAllard et"al] 
I2OOOI ). In addition, various photometric datasets are avail- 
able for the objects in Lupus. Optical and near-IR photom- 
etry in the Rc, Ic and Zwfi bands (at 0.6517, 0.7838 and 
0.9648 fim, respectively) was obtained from Comeron et al. 
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(2010, private communication). Tlie data were taken witli 
tiie Wide Field Imager (WFI) at tiie La Silla 2.2 m Tele- 
scope in Chile. Near-infrared Two Micron All- Sky Survey 
(2MASS) photometry at J, H and K (at 1.235, 1.662 and 
2.159 fim, respectively) is publicly available. Finally, Spitzer 
mid-infrared data at the IRAC and MIPS bands (at 3.6, 4.5, 
5. 8, 8.0, 24 and 70 /xm , respectively) are also available online 
Q l|Evans et al.ll2003h . 

For each object, its photometric data and correspond- 
ing NextGen model atmosphere are used (matching its 
spectral type and the gravitational acceleration for a pre- 
main-sequence star). The observed photometric data are 
corrected for extinction at all wavelengths from the ob- 
ject's visual extinction (Ay) , using the extinction law by 
IWeingartner fc Drainj (|200lh with Ry = 5.5. This extinc- 
tion is determined by SED fitting, with considering the de- 
rived spectral types. Subsequently, the NextGen atmosphere 
is scaled to the extinction corrected photometric data at a 
given band. The scaling band used is either the zwfi or J 
band, depending on the availability of the photometry and 
the quality of the fit. In Table [5J the extinction values and 
normalization bands are presented for all fitted objects. 

Figure [7] shows the constructed SEDs for 41 out of 54 
objects. Here, the solid black line is the NextGen model at- 
mosphere, the open squares are the observed photometry 
while filled circles show the extinction corrected data. SEDs 
could not be produced for 13 out of the 54 spectrally clas- 
sified objects. Either not enough infrared (both 2MASS as 
Spitzer) photometry is available, or the different photomet- 
ric data (optical, 2MASS and Spitzer) did not match well 
with each other. Eight of them are in crowded regions (18, 
19, 23, 26, 47, 48, 54, 89), which may have contaminated 
the photometry. None of the 13 objects are known variable 
stars. 

A wide variety of SEDs can be seen in Figure [T] Conse- 
quently, a wide variety in disk types is inferred. The infrared 
spectral index slope can be used as a characteristic to distin- 
guish between Class I, II and III objects, as defined by Lada. 
For the 39 objects with available data in K and MIPSl, the 
slope can be calculated as: 

«2-24nm = r T 

A24/^m A2/^m 

These values are shown in Table [S] One object has 
ct2-24)j.m consistent with a Class I source (q > 0) and 15 
consistent with Class III sources (a < —2), the majority 
of them are Class II, disk sources. Add itionally, object #6 7 
is a confirmed cold disk, as studied bv lMerfn et al.l l|2010l ). 
Cold disks are disks with large inner dust holes (i.e. lack of 
warm dust). Objects # 27, 79 and 80 have SEDs consistent 
with inner dust holes. Although confirmation of their na- 
ture as cold disks requires SED modelling that is out of the 
scope of this work, 4 out of 41 objects with SEDs in Fig- 
ure [7] amount to ~10% of cold disks in the sample, which is 
consistent with the range of 4 - 12% of cold disk frequency 
derived bv lMerm et^ l|2010l '). 
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5.1 Luminosities 

Stellar and disk luminosities can be derived for the objects 
for which SEDs could be plotted. Stellar luminosity is cal- 
culated with the formula: 

L. = 47rD^ J FxdX 

where D is the distance to the source and F\ is the stellar 
fiux - an integration of the NextGen model atmosphere 
scaled and normalized to the dereddened photometric data. 
Errors are derived from the uncertainty in the distance 
and flux, which includes an error in Ay of ±0.5 mag. The 
disk luminosity is the integrated excess emission above 
the stellar photosphere. Errors are derived from the stel- 
lar luminosity errors. Those two values are shown in TableO 



6 H-R DIAGRAM 

The age and mass of a star can be derived from theoreti- 
cal tracks, overlaid on physical H-R diagrams. The position 
of an object in the H-R diagram is determined by its tem- 
perature (derived in Section [4. 4p and luminosity (derived in 
Section 15. ip . 

6.1 Results 

The prc-main sequence evolutionary tracks of'Baraff e et al.l 
(1998, 2001) and Sicss ct al. (2000) are overlaid in the H-R 
diagrams where the objects in this sample are placed (Figure 
[8|. Thirteen objects are outside the range of the tracks (red 
and yellow circles in Figure [8l). 

Three objects (marked in red) are below the tracks: ID 
# 65, 71 and 78. By looking at their SEDs in Figure[7l it can 
be seen that those objects have substantial IR excess. The 
stellar fluxes are also not that high. These may be objects 
still surrounded by envelopes, which could be responsible for 
their low luminosity. Another explanation is that the objects 
are seen edge-on, such that the dust in the disk blocks some 
of the light of the star. In both cases, the calculated lumi- 
nosity will be too low, resulting in an incorrect placement 
in the HR diagram. Yet another reason may be that these 
are background sources and are in fact at a larger distance, 
which would increase the calculated luminosity. 

Ten objects (marked in yellow) are above the tracks: ID 
3, 13, 55, 60, 61, 63, 64, 68, 74 and 75. An easy explanation 
for the misplacement would be that the assumed distance is 
wrong and that the objects are actually closer to us. How- 
ever, even with very small distances (e.g. 60 pc), the objects 
are still above the tracks. Other effects have to play a role. 
None of them have detected lithium absorption and none 
of them, except object 13, show any infrared excess. One 
option is that they are part of a binary system that is un- 
resolved, resulting in a stellar luminosity that is too high. 
Another option is that they are evolved AGB-stars, which 
would explain the IR excess seen in object 13. The working 
assumption is that these objects, however, do not belong to 
the Lupus Clouds. 

For those 13 objects, no ages and masses were derived. 
For the other objects, the effective temperature, stellar and 
disk luminosity, age and mass can be found Table [3] As can 
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Figure 7. SEDs of all the objects. The black curve is the model atmosphere of a star at its given temperature. The open points are the 
observed optical, 2MASS and Spitzer data. The closed points are the extinction corrected data. Notice the wide variety in SEDs (and 
thus disk shapes). 



be seen from Figure [9l the two models agree on the range 
of ages and masses, albeit with a different distribution. The 
mean age is 3.6 Myr for the Baraffe models and 4.4 Myr 
for the Siess models. For the masses, both models deter- 
mine a mean stellar mass of onl y 0.3 M^. This res ult is in 
agreement with that derived bv lMerfn et al.l (|2008h . When 
deriving a luminosity function for the young stellar popula- 



tion in Lupus, the authors found it to peak at 0.2 Lq, which 
corresponds t o 0.2 Mq according to the PMS evolutionary 
tracks of .Baraffe et all lfi99& '). 
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Figure 8. HR diagram of the objects. Overlaid are the isochrones 
(blac k) and mass tracks fblue) . The top panel shows the tracks 
from lBarafFe et al.l lll998l . l200lh and the bottom panel the tracks 
from ISiess et alj feoOofT The corresponding ages (in Myr) are 
marked. The red and yellow circles are objects that do not follow 
the tracks. 



7 ACCRETION BASED ON Ha EMISSION 

From the data presented here, 45 of the objects (48 %) show 
the Ha hne in emission. The strength of the Ha emission hne 
is often used to distinguish between two classes of pre-MS 
stars: classical T Tauri stars (CTTS) and weak-line T Tauri 
stars (WTTS). CTTS are believed to correspond to stars 
that are actively accreting and WTTS are not, respectively. 
To this end, different methods have been proposed in the 
literature. Here, two methods are explored: the equivalent 
width and the full width of Ha at 10% of its peak intensity. 




5.0 5.5 6.0 6.5 7.0 7.5 -1.5 -1.0 -0.5 0.0 0.5 
og(Age/yr) log(M/MQ) 

Figure 9. Histograms of the ages and masses as derived from the 
two models. 



7.1 Ha Equivalent Width 

The value of equivalent width (EW) can be used to dis- 
tinguish between accreting and non-accreting objects, with 
the dividing threshold depe nding on the spectral type of the 
star. I White fc Basril (120031 ) proposed that a star is a CTTS 
if EW[Ha] ^ 3 A for K0-K5 stars, EW[Ha] ^ 10 A for K7- 
M2.5 stars, EW[Ha] ^ 20 A for M3-M5.5 stars and EW[Ha] 
> 40 A for M6-M7.5 stars. iBarrado v Navascues fc Martini 
(2003) extended this criterion to the substellar regime. Pos- 
itive values designate emission. It is important to note that 
stars with EW[Ha] values lower than these proposed lev- 
els are not necessarily WTTS. Extra confirmation is needed 
from other diagnostics such as the Lithium abundance. 

Following the criterion of White & Basri (2003), 24 of 
the 45 objects (53.3%) are considered to be actively accret- 
ing (values can be found in Table |4)|. For objects # 22, 49, 
50 and 59, no spectral type is available. However, since the 
EWs of these 4 objects are greater than 40 A, they are con- 
sidered accretors anyway. 



7.2 Full width of Ha at 10% of peak intensity 



IWhite fc Basril (|2003l ) proposed the use of another quantity 
to distinguish accreting and non-accreting objects: the full 
width of the line at 10% of the peak intensity, Ha[10%]. 
Lines with Ha [10%] > 270 km s~^ are, according to their 
measurements, assumed to be accretin g. Based on physical 
reason ing as well as empirical findings. Ijavawardhana et al.l 
l|2003h adopted 200 km s as a more reasonable accretion 
cut-off, especially for low-mass stars. This is the cut-off value 
used in this work. 

By performing a Gaussian fit (see Figure ITU)) . Ha [10%] 
values can be obtained. To account for the resolution, the 
measured FWHM of each profile must be deconvolved as- 
suming a Gaussian instrumental profile. For this sample, 
the Ha [10%] values can be found in Table U The errors 
are propagated from the error in the Gaussian fit to the 
er ror on the spectral re solution. According to the criterion 
of Ijavawardhana et al] (|2003), 25 of the 45 objects (55.5%) 
are considered accretors. The other 20 objects are not con- 
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Figure 10. Continuum subtracted, normalized profiles of the Ho emission lines (black line) overplotted with a Gaussian fit (red solid 
line). The red dotted line marks the centre of the Gaussian fit and the horizontal red line marks the width at 10% intensity. 
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Figure 11. Sample of continuum subtracted, normalized profiles 
of the Hq emission lines overplotted with a Gaussian fit (red 
dashed line) and a Voigt fit (thick blue solid line). The red dotted 
line marks the center of the Gaussian fit. 
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Figure 12. Histogram of the mass accretion rates derived from 
the Ha [10%] value. The shad ed area are o bjects out of the range 
for which the relationship of iNatta et al.l (2004) was calibrated. 
Those objects are generally classified as non-accretors. 



sistent with CTTS, even within their uncertainties. They are 
classified as non-accretors. 

Some line profiles show broad wings at the bottom, 
which cannot be properly fitted by a Gaussian. In order 
to reproduce those lines, Voigt profiles were fitted to all the 
lines. In most cases, the fit is either the same or slightly bet- 
ter at the wings as can be see in Figure 1111 However, this 
correction has no significant effect on the full width at 10%. 
Due to this insignificant difference, only the values from the 
Gaussian fits are further considered in this work. 

The EW and Ha[10%] methods are generally consistent 
in the classification of sources. Only for three objects (#65, 
80, 81) do the two methods give discrepant results. This may 
be the result of the different cut-offs in the EW method, 
depending on the spectral type. 



7.3 Mass accretion rate 



The Ha[10%] value can not only be used as an indicator 
of accretion, but al so allows a quantit ative estimate of the 
mass accretion rate. lNatta et al] l|2004 ) derived a correlation 
between Ha [10%] and mass accretion rate Mac (derived from 
veiling and Hq profile model fittings): 

\og Mac = -12.89(±0.3) + 9.7(±0.7) • 10~^Ha[10%] 

where Mac is in Mr yr~^ and Hq[1 0%] in km s"^ The 
relation derived by iNatta et al.l (|2004l ) was only calibrated 
for -11 < log A/ < -6. The cut-off of 200 km s"^ agrees 
with a cut-off of lO-"-^'^ M© yr"^ for M. With this relation, 
mass accretion rates can be calculated for the objects in 
this sample with Ha in emission (see Table Q. In Figure 
1121 the distribution of the mass accretion rates is shown. 
Typical uncertainties are an order of magnitude. Most of the 
accretors have mass accretion rates that are consis tent with 
the typical range for classical T Tauri stars (Na tta et al.l 
12004 ISiciHa-Aguilar et a"i]|2010l ). 

A fit of the mass accretion rate as a function of the 
stellar mass is attempted. This fit reveals that a Mac oc 
relation is valid for both models, within the large errors. 



The same relation, and huge scatter, can be f ound in the 
literature (|Natta et al.ll2005l : lATcala et al.ll201ll ). 



8 DISCUSSION 

Our work can be added to some of the objects classified 
as young stars bv lMerm et al.l (|2008l ) that have previously 
bee n spectrally classified in t he literature in different stud- 
ies llL6pez Marti et ai]|2005l: lAllers et a.l.|[200^ : lAllen et al] 
I2OO7I : iMerfn et al.ll2007l : IComer6nll2008l ). This amounts to 
91 of t he 159 IR excess sources in the sample of lMerfn et al.l 
(|2008l) to have known spectral types. This fairly complete 
sample can be compared to other nearby star-forming re- 
gion s with equally c ompl ete stellar characterization. 

iLuhman et"al] (2003) studied the young cluster IC348. 
Spectral types are determined with lo w-resolut i on op tical 
spectroscopy for 169 objects. Similarly, iLuhmanI (|2007l ) de- 
termined spectral types for the 85 targets in the star form- 
ing region Chamaeleon I. Figure [TS] shows the distribution 
of spectral types in Lupus (this work + literature), in IC 
348 and Chamaeleon I. The vertical dashed line shows the 
completeness of the IC 348 and C ha I samples , which is com- 
parable to the limit achieved bv lMerfn et al.l l|2008l ). It can 
be seen that the three regions show very similar distribu- 
tions of stellar types, from which is inferred that the IMF of 

Lupus is very similar to those of IC 348 and C ha I. 

It is important to note that the sample of lMerfn et al.l 

(|2008| ) is selected on the basis of IR excess and it is, by 
definition, biased against young stars without disks (i.e. no 
IR excess). A similar study characterizing the young stellar 
population without disks in Lupus is still lacking. 

When comparing the stellar luminosity with the disk 
luminosity, a correlation Ldisk oc L* can be found with a 
correlation coefficient of 0.82. The spectral index a2-24/jm 
is also compared with the stellar parameters (mass, age, lu- 
minosity). Slight trends can be seen in that the IR excess 
increases with mass and decreases with the stellar and disk 
luminosity. However, the scatter is huge in all the fits. More 
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Figure 13. Distribution of spectral types in Lupus (black solid 
line), IC 348 (green dot-dashed line) and Chamaeleon I (blue 
dashed line). For Lupus, the distribution includes both our results 
as well as literature spectral types for sources not in our sample. 



observations and a larger sample are needed to make signif- 
icant conclusions. 



9 CONCLUSIONS 

We have presented a spectroscopic survey at optical wave- 
lengths designed to determine spectral types and confirm 
the pre-main sequence nature of a sample of young stellar 
objects in the Lupus Clouds found in the Spitzer c2d survey. 

• Spectral types were determined for 54 stars belonging 
to Lupus I, III or IV. The sample consists mostly of M- 
type stars (90%), but also a few K-type stars. No early type 
object was found in this sample. The distribution of spectral 
types peaks at M4-M6. The distribution of spectral types is 
very similar to that of well studied star-forming regions like 
IC 348 and Chamaeleon I. 

• The objects were placed in a H-R diagram after effec- 
tive temperatures and luminosities were derived. Compari- 
son with theoretica l isoch rone s and mass t racks from mod- 
els of iBaraffe eFaP (l998') and lSiess etlo] l|2000l ) yield indi- 
vidual ages and masses for the objects. 10 objects are too 
luminous to belong to the clouds, while 3 objects are too 
faint. The very faint objects could still belong to the clouds 
by having a remnant envelope or an edge-on disk dimming 
some of the stellar luminosity. 

• The theoret ical models by iBaraffe et al.l (|l998l ) and 
ISiess et al] (|2000l ) imply a population of YSOs concentrated 
in the age range between 1 and 5 Myr. The mean age is f ound 
to be 3.6 and 4.4 Myr with the IBaraffe et al] l|l998h and 
ISiess et al.l l|2000l ) tracks, respectively. Individual masses 
range from 0.1 to 1.0 Mq, with mean values of 0.3 Mq for 
both models. 

• About half of the sample shows the Ha line in emission, 
an important indicator of accretion. This relationship was 
explored in two different ways: the equivalent width of Ha, 
and its full width at 10% of peak intensity. This confirms 25 
objects (or 56% of the YSO sample) to be actively accreting 
objects classified as classical T Tauri stars. The quantitative 




estimate of the mass accretion rate Mac based on the full 
width of Ha at 10% of the peak intensity yields a broad 
distribution of values (~ 10~^^ — 10~* M0yr~^), typical of 
T Tauri stars. 
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Table 1: ID, name, cloud, position, observing dates and exposure times 
of the objects. 



ID 


Name 


Cloud 


RA 
(des 


?) 


Dec 
(deg 


) 


Observing Date 


(s) 


Mark^ 


1 


2MASSJ16075475-3915446 


Lupus III 


241 


9780 


-39 


2623 


2008-02- 


21 


630.0 


u 


2 


2MASSJ16080175-3912316 


Lupus III 


242 


0070 


-39 


2087 


2008-02- 


21 


630.0 


u 


3 


2MASSJ16080618-3912225 


Lupus III 


242 


0260 


-39 


2062 


2008-02- 


21 


630.0 


G 


4 


2MASSJ16081497-3857145 


Lupus III 


242 


0620 


-38 


9541 


2008-02- 


21 


630.0 


u 


5 


2MASSJ16084747-3905087 


Lupus III 


242 


1980 


-39 


0856 


2008-02- 


21 


630.0 


u 


6 


2MASSJ16085373-3914367 


Lupus III 


242 


2240 


-39 


2436 


2008-02- 


21 


630.0 


u 


7 


2MASSJ16085529-3848481 


Lupus III 


242 


2300 


-38 


8134 


2008-02- 


22 


1000.0 


G 


8 


2MASSJ16085953-3856275 


Lupus III 


242 


2480 


-38 


9411 


2008-02- 


21 


630.0 


u 


9 


AKC2006-17 


Lupus I 


234 


8640 


-34 


8121 


2008-02- 


23 


800.0 


u 


10 


IRACJ16083010-3922592 


Lupus III 


242 


1250 


-39 


3831 


2008-02- 


26 


1000.0 


u 


11 


IRACJ16084679-3902074 


Lupus III 


242 


1950 


-39 


0354 


2008-02- 


21 


630.0 


u 


12 


IRACJ16093418-3915127 


Lupus III 


242 


3920 


-39 


2535 


2008-02- 


26 


1000.0 


u 


13 


IRAS15567-4141 


Lupus IV 


240 


0310 


-41 


8302 


2008-02- 


21 


550.0 


G 


14 


IRAS15589-4132 


Lupus IV 


240 


5900 


-41 


6816 


2008-02- 


24 


1000.0 





16 


Lup604s 


Lupus III 


242 


0010 


-39 


0497 


2008-02- 


21 


630.0 


G 


17 


Lup607 


Lupus III 


242 


1170 


-39 


2194 


2008-02- 


21 


630.0 


G 


18 


Lup608s 


Lupus III 


242 


2850 


-39 


0619 


2008-02- 


21 


630.0 


G 


19 


Lup617 


Lupus III 


242 


2010 


—39 


1554 


2008-02- 


21 


630.0 


G 


20 


Lup650 


Lupus III 


242 


4580 


—38 


8174 


2008-02- 


22 


1000.0 


G 


21 


Lup654 


Lupus III 


241 


8480 


-39 


0861 


2008-02- 


21 


630.0 


F 


22 


Lup706 


Lupus III 


242 


1560 


—39 


3864 


2008-02- 


26 


1000.0 


F 


23 


Lup710 


Lupus III 


242 


3210 


—39 


4527 


2008-02- 


23 


1000.0 


G 


24 


Lup714 


Lupus III 


241 


9950 


—39 


4097 


2008-02- 


26 


1000.0 


G 


25 


Lup802s 


Lupus III 


242 


9630 


-38 


8514 


2008-02- 


26 


1000.0 


G 


26 


LnnSlOs 


Lupus III 


242 


4770 


—39 


2009 


2008-02- 


22 


1000.0 


G 


27 


LupSlSs 


Lupus III 


242 


4850 


-38 


9975 


2008-02- 


22 


1000.0 


G 


28 


LupusSMMS 


Lupus III 


242 


3250 


-39 


0815 


2008-02- 


24 


10.0 


F 
















2008-02- 


26 


1000.0 




29 


NTO2000-0526. 9-5630 


Lupus III 


242 


1170 


-39 


0735 


2008-02- 


21 


630.0 


u 


30 


NTO2000-0532. 1-5616 


Lupus III 


242 


2230 


-39 


0692 


2008-02- 


21 


630.0 


u 


31 


NTO2000-0536. 7-5943 


Lupus III 


242 


2430 


-39 


1265 


2008-02- 


24 


1000.0 


F 


32 


NTO2000-0536. 7-5956 


Lupus III 


242 


2430 


-39 


1304 


2008-02- 


24 


10.0 


F 
















2008-02- 


26 


1000.0 




33 


NTO2000-0537.4-5653 


Lupus III 


242 


2460 


-39 


0794 


2008-02- 


21 


630.0 


u 


34 


NTO2000-0540.9-5757 


Lupus III 


242 


2600 


-39 


0970 


2008-02- 


21 


630.0 


F 


35 


NTO2000-0546. 4-5934 


Lupus III 


242 


2830 


—39 


1241 


2008-02- 


24 


10.0 


F 
















2008-02- 


26 


1000.0 




36 


NTO2000-0554.9-5651 


Lupus III 


242 


3180 


-39 


0788 


2008-02- 


24 


1000.0 


F 


37 


NTO2000-0558.8-5610 


Lupus III 


242 


3350 


-39 


0671 


2008-02- 


24 


10.0 


F 
















2008-02- 


26 


1000.0 




38 


NTO2000-0601. 7-5616 


Lupus III 


242 


3460 


-39 


0687 


2008-02- 


24 


10.0 


F 
















2008-02- 


26 


1000.0 




39 


NTO2000-0605.1-5606 


Lupus III 


242 


3610 


-39 


0660 


2008-02- 


26 


1000.0 


U 


40 


NTO2000-0605.6-5437 


Lupus III 


242 


3630 


-39 


0412 


2008-02- 


21 


630.0 


U 


41 


NTO2000-0614.0-5414 


Lupus III 


242 


3980 


-39 


0349 


2008-02- 


21 


630.0 


u 


42 


NTO2000-0615.6-5616 


Lupus III 


242 


4050 


-39 


0686 


2008-02- 


24 


1000.0 


u 


43 


NTO2000-0615.6-5953 


Lupus III 


242 


4050 


-39 


1291 


2008-02- 


24 


1000.0 


u 


44 


NTO2000-0615. 8-5734 


Lupus III 


242 


4060 


-39 


0905 


2008-02- 


24 


1000.0 


u 


45 


NTO2000-0617.7-5641 


Lupus III 


242 


4140 


-39 


0755 


2008-02- 


24 


1000.0 


u 


46 


NTO2000-0619.6-5414 


Lupus III 


242 


4210 


-39 


0349 


2008-02- 


21 


630.0 


u 


47 


Par-Lup3-1 


Lupus III 


242 


0670 


-39 


0511 


2008-02- 


21 


630.0 


G 


48 


Par-Lup3-2 


Lupus III 


242 


1490 


-39 


0632 


2008-02- 


21 


630.0 


G 


49 


Par-Lup3-3 


Lupus III 


242 


2060 


-39 


0942 


2008-02- 


21 


630.0 


F 


50 


Par-Lup3-4 


Lupus III 


242 


2140 


-39 


0917 


2008-02- 


21 


630.0 


F 
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Table 1 - Continued 



ID 


Name 


Cloud 


RA 




Dec 




Observing Date 


texp 


Mark^ 










(deg) 


(deg) 




(s) 




01 


qC'p T iiri^ 1 
DO ±-ijUpo-l 


Lupus 


III 


242, 


.9990 


—38, 


.3940 


ZUUo-UZ-ZO 


1 000 n 

lUUU.U 






oo ± (jZQj iOoyzo.z-4:Zoou i 


Lupus 


IV 


239, 


.8550 


-42, 


.5853 


9nns 09 9fi 


1 000 

lUUU.U 


a 
\j 


Oo 


oo J- CZCIJ iOUUUU.O-^ZZlOo 


Lupus 


IV 


240 


.0020 


—42, 


.3659 


9nnR 09 9f^ 


1 000 n 

lUUU.U 


yj 


04 


QQTi^9rl 71 (^nnn9 /I /19991fi 
oo 1 CZQJ iDUUUZ.4-4ZZZ10 


Lupus 


IV 


240, 


.0100 


—42, 


.3708 


ofins n9 9f^ 


1 nnn o 

iUUU.U 


n 
yj 


OO 


'^^Tr9H 71 f^O1 1 1 zLI ^7^0 


Lupus 


IV 


240, 


.2980 


-41 


.6250 


onnc 09 91 


oou.u 


n 


oo 


OO J- CZQJ lOUl'yiO.o-^loDUD 


Lupus 


IV 


240, 


.4300 


—41, 


.6016 


900ft 09 91 


oou.u 


n 


1 


QCITr-Orl 71 fin99Q Q A1 t^l 1 1 
OO iczujiouzzy.y-^ioiii 


Lupus 


IV 


240 


.6250 


—41 


.8530 


9008 09 94 


1 000 
lUUU.U 


p 


oo 


^'^Tr9H 71 ^070'^ Q "^Ql 119 
oo ± L-zLij iuu / uo.y-oyiiiz 


Lupus 


III 


241 


.7660 


-39 


.1865 


9008 09 91 


(S'^O 
uou.u 


n 
\j 


oy 


oo 1 CZQJ lou ( uo.o-oy ( 


Lupus 


III 


241 


.7860 


—39 


.2354 


900S 09 91 


fi^O 
OoU.U 


r 


fin 


oo ± czaj ±ou / oo.o-oyu / io 


Lupus 


III 


241, 


.9800 


—39, 


.1215 


9008 09 91 


fiQO 
OoU.U 




U-L 


OO ± CZU.J xuuouo.u~ooozzy 


Lupus 


III 


242, 


.0130 


—38, 


.8749 


900S 09 94 
zuuo~uz-z^ 


1 000 

lUUU.U 






QC;Tr>9rl 71 fiflQni A ^Q9f^19 

oo ±czajiouyui.4-oyzoiz 


Lupus 


III 


242, 


.2560 


—39, 


.4200 


900S 09 9^ 


1 000 
lUUU.U 


n 

or 


Oo 


'^c;T'r>9rl 71 fiflQ'^/l 1 '^QI'^/19 


Lupus 


III 


242, 


.3920 


-39, 


.2284 


900R 09 91 


DOU.U 


\j 




QciTr^9ri 71 fil nnn i '^sf^j^ni 

oo ± CZQJ lOlUUU. i-ooO-iUi 


Lupus 


III 


242, 


.5000 


-38, 


.9001 


900ft 09 99 
zuuo-uz-zz 


1 000 
lUUU.U 


r< 

or 


OO 


QC;'Po9rl 71 fil ni ^ 1 ^8Afi1 7 
oo J- CZQJ lOlUlO. 1-0041:01 f 


Lupus 


III 


242, 


.5540 


—38, 


.7712 


900S 09 99 
zuuo-uz-zz 


1 000 
lUUU.U 


n 

Vjr 


66 


SSTc2dJ161019.8-383607 


Lupus 


III 


242, 


.5830 


—38, 


.6018 


2008-02-23 


1000.0 


G 


67 


SSTc2dJ161029. 6-392215 


Lupus 


III 


242 


.6230 


-39, 


.3707 


2008-02-23 


1000.0 


C 


Oo 


oo 1 CZQJ lOlUo'i.O-OOl^OU 


Lupus 


III 


242 


.6440 


—38, 


.2473 


9008 09 9f^ 
ZUUo-UZ-ZO 


1 000 n 

lUUU.U 


or 


oy 


QC;Td^9rl 71 fil 1 ^1 Q "^811 in 

oo 1 CZQJ loiioi.y-ooiiiu 


Lupus 


III 


242, 


.8830 


—38, 


.1861 


9008 09 9 
ZUUO-UZ-ZO 


1 000 

lUUU.U 


r< 

or 


70 


ClClTn9rl 71fi1 1AA Q ^8^94 
OO 1 CZQJ ±oi i^^.y-oooz^o 


Lupus 


III 


242, 


.9370 


—38, 


.5458 


9008 09 9^ 
zuuo uz-.^o 


1 000 

lUUU.U 


Or 


71 


^^Tr9H 71 f^l 1 zl8 7 "^Rl 

oo 1 CZLIJ ±U114:0- i -Ool i OO 


Lupus 


III 


242, 


.9530 


-38, 


.2994 


9008 09 9^ 

ZUUO-UZ-ZO 


1 000 n 

lUUU-U 


n 

Or 


79 


c;QT'r'9r1 71 fil 90/1 ^ '^SOQf^Q 

OO 1 CZQJ loizu^.o-oouyoy 


Lupus 


III 


243, 


.0190 


-38, 


.1663 


9008 09 9f^ 
ZUUO-UZ-ZO 


1 000 n 

lUUU.U 


TT 

U 


7*^ 


^'^Tr9H 71 fil 91 1 9 '^8'H990 

oo 1 CZLIJ ±U-LZ±±.Z-OOOZZU 


Lupus 


III 


243, 


.0470 


—38, 


.5389 


9008 09 99 
zuuo-uz-zz 


1 000 n 

lUUU-U 


TT 
U 


74 


SSTc2dJ161219.6-383742 


Lupus 


III 


243, 


.0820 


-38, 


.6283 


2008-02-22 


1000.0 


G 


75 


SSTc2dJ161251.7-384216 


Lupus 


III 


243, 


.2150 


-38, 


.7044 


2008-02-22 


1000.0 


G 


/ u 


oz±uu 


Lupus 


III 


242, 


.1070 


—39, 


.1003 


9008 09 91 
zuuo-uz-z± 

2008-02-23 
2008-02-26 


■^00 
oou.u 

200.0 

200.0 


Or 


77 




Lupus 


III 


242, 


.1180 


-39, 


.0922 


9008 09 91 
zuuo-uz-z ± 


fi^O 
uou.u 


Or 


7S 




Lupus 


III 


242, 


.1240 


—39, 


.0530 


9008 09 91 

ZUUO UZ-Z L 


R^O 
uou.u 


Or 


79 


Szl03 


Lupus 


III 


242 


.1260 


—39, 


.1030 


2008-02-21 


630.0 


G 


80 


Szl04 


Lupus 


III 


242, 


.1280 


—39, 


.0968 


2008-02-21 


630.0 


G 


81 


<^'7^ 07 

OZIU 1 


Lupus 


III 


242, 


.1740 


—39, 


.0269 


9008 09 91 
ZUUO-UZ-Z 1 


fi^O 
OoU.U 


r< 

or 


89 


Qr7l 087^ 

OZIUOO 


Lupus 


III 


242, 


.1790 


—39, 


.1040 


9008 09 91 
ZUUO-UZ-Z 1 

9008 09 9^ 
ZUUO uz-zo 

9008 09 ^f, 
zuuo-uz-zu 


fi^O 
OoU.U 

900 
zuu.u 

900 
zuu.u 


P 
or 


oO 


ci-T-i 1 n 

OZllU 


Lupus 


III 


242, 


.2150 


-39, 


.0549 


9008 09 91 
ZUUO-UZ-Z 1 

9008 09 9^ 
ZUUO uz-zo 

2008-02-26 


•^00 
oUU.U 

900 
zuu.u 

200.0 


P 
or 


84 


Szll3 


Lupus 


III 


242, 


.2410 


-39, 


.0396 


2008-02-21 


630.0 


G 


8^; 

oO 


Cl-yl 1 A 

OZll'i 


Lupus 


III 


242, 


.2580 


—39, 


.0867 


9008 09 91 
ZUUO-UZ-Z 1 

2008-02-23 
2008-02-26 


•^00 
oUU.U 

200.0 

200.0 


p 

or 


8fi 
oO 


Qrzl 1 8 

OZl lo 


Lupus 


III 


242, 


.4530 


-39, 


.1880 


9008 09 99 
ZUUO-UZ-ZZ 


1 000 
lUUU.U 


P 
or 


87 


SzlSO 


Lupus 


IV 


240, 


.1290 


-41, 


.7269 


2008-02-21 


550.0 


G 


OO 
OO 




Lupus 


IV 


240, 


.8730 


-41, 


.6673 


onno no 0/i 
zUUo-Uz-z4 


1 nnn n 
lUUU.O 


/-I 
(j 


89 


Sz94 


Lupus 


III 


241, 


.9570 


-39, 


.0746 


2008-02-21 


630.0 


G 


90 


Sz96 


Lupus 


III 


242, 


.0520 


-39, 


.1425 


2008-02-21 


630.0 


G 


91 


Sz97 


Lupus 


III 


242, 


.0910 


-39, 


.0726 


2008-02-21 
2008-02-23 
2008-02-26 


630.0 
200.0 
200.0 


G 


92 


Sz98 


Lupus 


III 


242, 


.0940 


-39, 


.0794 


2008-02-21 
2008-02-23 

2008-02-26 


300.0 
200.0 

200.0 


G 


93 


Sz99 


Lupus 


III 


242, 


.1000 


-39, 


.0970 


2008-02-21 
2008-02-23 


630.0 
200.0 


G 



A. Mortier et al. 



Table 1 - Continued 



ID 


Name 


Cloud 


RA 


Dec 


Observing Date 


tcxp Mark^ 








(deg) 


(deg) 




(s) 












2008-02-26 


200.0 



Table 2: Spectral types, visual extinction (Av), normalization bands and 
the spectral index, calculated between the 2MASS K-band (2 ^.m) and 
the MIPSl band (24 /im), for the classifiable objects. Objects for which 
SEDs could not be produced do not have all values. The last column 
indicates if Lithium was detected. 



ID 


SpT 


Range 


Av 


Band 


Q!2-24fim 


Lithium? 


3 


M6.5 


M4 - MS. 5 


5.30 


J 


-1.91 




7 


M3 


M2 - M4 


2.50 


Z 


-1.25 




13 


M6.5 


M4 - MS. 5 


2.00 


J 


-1.54 




16 


M5.5 


M4 - M6 


1.20 


Z 


-1.21 


yes 


17 


M5.5 


M4 - M6 


0.00 


Z 






18 


M5.5 


M4 - MS. 5 










19 


M6 


M4 - MS. 5 










20 


M3.5 


M2 - M4 










23 


M4.5 


M3 - M6 








yes 


24 


M6 


M4 - MS. 5 








yes 


25 


M5.5 


M4 - M6 










26 


KO 


KO - K2 










27 


M6 


M4 - MS. 5 


1.40 


Z 


-1.01 




47 
48 


M5.5 
M5 


M4 - M6 
M4 - M6 










51 


M6 


M4 - MS. 5 


1.80 


J 


-1.08 


yes 


52 


M4.5 


M3 - M6 








yes 


54 


M3.5 


M2 - M4 








yes 


55 


MS. 5 


M6 - M9.5 


0.00 


J 


-2.30 




56 


M6 


M4 - MS. 5 


1.20 


J 


-2.36 




57 


M7 


M4 - M9 


2.60 


J 


-2.27 




58 


M5.5 


M3 - M6 


2.20 


Z 


-0.17 




60 


M9 


MS. 5 - M9.5 


S.OO 


Z 


-2.16 




61 


M5.5 


M4 - M6 


2.00 


Z 


-2.39 




62 


M3.5 


M2 - M4 


0.60 


J 


-1.10 


yes 


63 


M7.5 


M6 - M9 


3.40 


Z 


-2.49 




64 


M6.5 


M6 - MS. 5 


2.80 


J 


-1.82 




65 


M4 


M3 - M6 


0.70 


Z 


-0.21 




66 


M6 


M4 - MS. 5 


0.40 


J 


-1.18 


yes 


67 


M5.5 


M3 - M6 


0.00 


J 


-0.93 


yes 


68 


M6.5 


M6 - MS. 5 


2.90 


J 


-2.31 




69 


M2.5 


MO. 5 - M4 


1.90 


J 


-2.30 




70 


M6.5 


M4 - MS. 5 


0.00 


J 


-1.99 


yes 


71 


K2 


KO - K7 


1.60 


J 


-0.60 




74 


MS 


M6 - MS. 5 


1.60 


J 


-2.34 




75 


MS. 5 


M6 - M9.5 


1.00 


J 


-2.15 




76 


M4.5 


M3 - M6 


1.20 


Z 


-0.80 


yes 


77 


M4.5 


M3 - M6 


0.70 


Z 


-1.68 


yes 


78 


K2 


KO - K7 


2.50 


J 


0.59 




Continued on next page 



means good spectra, U non-detections, F featureless spectra, and O not useful spectra. See 14.31 for more details. 



YSOs in the Lupus Molecular Clouds 



Tabic 2 Continued 





ID 


SpT 


Range 


Av Band 
(mag) 




Lithium? 






7Q 


M4.5 


M3 - M6 


1 nn 7 

l.UU Zj 


— U.oO 


yes 






SO 


M5.5 


M3 - M6 


n 70 7 




yes 






ol 


M6.5 M4 - M8.5 


n on 7 










oz 


M5.5 


M3 - M6 


n sn 7 

U.OU Zj 




yes 






83 


M3 M0.5 - M4 


0.20 J 


-0.71 


yes 






84 


M1.5 


K7- M2 


2.40 Z 


—0.67 








85 


M4 


M3 - M6 


1.30 Z 


-0.62 


yes 






86 


K7 


K2 - MO 


2.60 J 


-0.90 


yes 






87 


M2 MO. 5 - M4 


0.60 J 


-0.94 


yes 






88 


K2 


KO- K7 






yes 






89 


M4.5 


M3 - M6 












90 


M2 


MO - M3 


1.43 J 


-0.91 


yes 






m 
yi 


M4.5 


M3 - M6 


0.30 Z 


-1.24 


yes 






Q9 


K5 


K2 - MO. 5 


2.50 J 


-0.64 


yes 






93 


M4 


M3 - M6 


0.40 J 


-1.12 


yes 






Table 3: 


: Effective temperature Tefi, stellar and disk luminosity, L« and 






Ldisk, 


;e and mass for Baraffe et al. (1998 - 2001) and Siess et al. (2000). 




ID 


Teff 


L, 


Ldisk 


Ages 


Masss 


Ages 


Masss 




(K) 


(Lo) 


(Le) 


(Myr) 


(Mo) 


(Myr) 




3 


OQQC: r|+336.0 


9 nq+°-^° 


0.610l°;« 










7 


3415.0tl«° 


06+° °3 
'^•'-'"-0.02 


o.oiol°:°l 


13.64li*2f 


n 0C+O.I6 
U.OO_o 12 


9-901^:!^ 


27+0 07 
"•^'-o.oe 


13 


2935.01335.0 


^- '^-0.94 


1.930lJ:?^ 










16 


3057.5le7Y' 


"•'-'^-0.03 


o.o2oi°:;g 


1 -,7+2.79 

-0.78 


l2+°-i3 


3.85l°3:°^ 


1 4+0.08 
"•1^-0.02 


17 


3057.5lg7;^/' 


m+°°° 

^•'-'1-0.00 


o.oool°:»j 


12.75lfif 


0.091°:°° 


16.25l*«4f 


o.o9l°:««2 


27 


2990.01435.0 


n 04+°°^ 

U-U'±_o.02 


0.0101°:- 


1 07+7. 16 
-'-•0'-1.87 


0.09l°:J3 


^•^■5-2.82 


o.o9l°:J° 


51 


2990.0143";!; 


12+°°* 

^•1^-0.03 


0.0101°:- 




o.iil°:i? 


ft2+0-50 
•J.uz_0 59 


n iy+0.06 
'J-i'l -0.06 


55 


2555.0l«-; 


1 7+0-08 


o.o5ol°:°| 










56 


2990.0l^»»;0 


S7+°-i^ 


o.oGoi;;:;;^ 


00+° " 

U.UU_o 00 


o.i8l°:i^ 


^•S-L-LSl 


"•^"-0.26 


57 


2880.011°° 


1 S+°°^ 


o.o3oi°:°l 


00+2 ** 
"J-'^'^-o.oo 


06+5 02 
'-'•'^"-0.00 


4n+-3.40 
'J-'±U_3 40 


n 24+-0-24 

'J-^*-0.24 


58 


3057.51^^^5'^ 


01+° °^ 


o.o2ol°:°i 


9.831^:^°^ 


09+°-^° 

U.Uy_o 02 




09+" "' 

"•O»_o.o2 


DU 


2400.0Ii(j(j 


1 Q4+0-95 


o.8ool°:!J 










61 


3057.5167^5^ 


1 99+0.60 
J-.^^_0.42 


o.i9ol°:?3 










62 


3342.5l72''5^ 


16+006 


o.o6ol°:°3 


2-86lg:« 


SI +0-20 


4Q+1.67 
^■^^-0.67 


r. 97+0.09 

"•^'-o.oa 


63 


2795.0l^«^:° 


9 qo+1.43 


i.i4ol°:«t 










64 


2935.OI38Q g 


1 rc+0.60 
1.0D_o 45 


0.4501°:^* 










65 


3270.0l^fo:g 


0.001°:°° 


o.oool°:°i 










66 


2990.01^*°:° 


o.o3l°:°l 


o.oool°:°J 


3.07li^of 


09+0- 12 

U.U»_o 09 


O.UZ_2 9]^ 


o.o8l°:J° 



continued on next page 



without ages and masses do not follow the evolutionary tracks. 



Mortier et al. 



ID Teff L« Ldisk Ages Masss Ages Masss 





(K) 




(L0) 




(L©) 


(Myr) 




(Mo) 


(Myr) 




(Mo) 


67 


3057.5lg"5^ 


0, 


QQ-l-u.uy: 
■•^^-0.03 


0, 


.0201°:°? 


i-i7lo:?8 


0. 


1 + 0. 24 
■^^-0.02 


^„_|_0 Of) 

O.OJ_3 g5 


0, 


-1 A -1-0 1 4 
■^^-0.02 


68 


2935.0l3gQ g 


7, 


■1^-2.07 


0, 


•990la82 














69 


3487.51^^1^:^ 


0. 


-1 Q-l-0.05 
•-'-'^-0.04 


0, 


^i-.rt-1-O 02 

.0201°:°^ 




0. 


A A-\-0 26 
■^^-0.18 


_ rr J -1-8 90 

5.741?:?° 


0. 


„„_|_0.19 
•<J<J_0.09 


70 


2935.0l^3«:° 


0. 


j-i/i-l-O 02 


0, 


^n^-l-O 01 

■oool°:°j 


r\ A A-\-7 07 

0.44lJ:°I 


0. 


*-in-l-0 15 

■UO-0.08 


-2.32 


0. 


r»n-l-0 12 

•"°-0.08 


71 


4900.0l|°°:° 


0. 


r\ A-\-0 02 


0, 


,-.r».-.-l-0 01 

.0201°:°; 














74 


2710.01?=°:° 


2. 


A A-\-0 93 

.441°:?? 


0, 


■380l°:^| 














75 


ZDDD. U_255.o 


2. 


•3olo:i7 


0, 


f-i-T^-l-0 47 














76 


3197.5^207.5 


0. 


0^7-1-0.13 
■^'-0.09 


0, 


■lool^:^^ 


n co+1-73 
U.DO_o_53 


0. 


i-.^-1-O 20 


-\ r\ A -1-0 75 

1.941?:^° 


0. 


90+0.08 
•^''-0.02 


77 


--1-21 7 

3197.5l^j;:^ 


0, 


, ^ _i_n on 
■^^-0.14 


0, 


•06ol°:J° 


„ ^ --I-1 25 


0. 


■^ ^-1-0 2fi 

■-■-^-0.01 


U-O^_0'.82 


0. 


9c:+0.07 
•^^-0.02 


78 


4900.0lg°:° 


0, 


/-lO-l-O 01 


0, 


1 or»-l-0 04 

•1301°:°;^ 














79 


3197.5l^J^:^ 


0, 


r»r\+{) 10 





.0601;;:;;;^ 


r\ r\r\-\-l 86 

QQ ' ^ 


0, 


nr» + 18 

■2oi;;:j? 


. n + 67 

2.481°:;;^ 


0, 


99+0.09 
■^^-0.04 


80 


3057.515^5^ 


0, 


^ Q -1-0.06 
■-'-'^-0.04 





.0401;;:;;^ 


„ , . -1-5 03 

0-641;;:;;? 


0, 


-| Q + 0.24 
■-L-^-O-Ol 


f-* r-\-^-l-1 Q6 


0, 


ic+O.M 
■l'-'-0.02 


81 


2935.01^^^.0 


0, 


H 1- -l-ft 07 

■-•-^-o.os 





.oioi;;:;;t 


o-oolf,:gg 


0, 


■-'-'-'-0.10 


Q 00 H 3. 28 
3-2»_3,28 




82 


3057.515";"' 


0, 







.0301°:;;^ 




0, 


1 + 0. 24 
■-'-^-0.02 


3-861^;;;; 


0, 


1 4+0.14 
■1^-0.02 


83 


3415.01?^:;; 


0, 


■^"-0.08 


0, 


.0901°:°^ 


n F--1-7 10 

^•<=>^_1.20 


0. 


. f^-l-0 38 

■^'-'-0.12 


^- ' "^-0.38 


0, 


qi +0.20 

■•Ji-o.oe 


84 


3632.5172^5^^ 


0, 


+0.05 
■-'^J^-0.04 


0, 


rv A rv-l-0 03 

.0401°:°^ 


-1 <-i ^ -1-59 62 

19.64l5°g6=^ 


0. 


rr>-l-0 08 

■^°-0.08 


n nr+5() 10 

9-95l^'9^° 


0, 


4Q+0.24 
■^'-'-0.04 


85 


3270.0l^«:° 


0, 


40+0-24 


0, 


.2201°:;? 


'J-^D-0.26 


0. 


Q9+0. 16 
■"^^-0.10 


1 19+0.96 
J^-J^^-1.12 


0, 


97+0.05 
-0.01 


86 


4060.0ll«:° 


0, 


09+0.35 
■^^-0.27 


0, 


.5601°:- 


.1 79+20.71 
■J- ' ^-1.37 


1. 


1 e+0.()4 
■J^"-0.16 


9 c:4+13.53 


0, 


■ '"-0.20 


87 


3560.01^1^:^ 


0. 


91+0.10 
•^-^-0.07 


0, 


.060l°:°« 


cr 77+7.30 
^- ' ' -4.20 


0. 


C9+0.23 
■^^-0.25 


3.741?:°^ 


0. 


07+0.14 
••^'-0.12 


90 


3560.0l?«°:° 


0. 


09+0.31 
•°^-0.24 


0, 


.27oi°:;« 


0-88lJ:^« 


0. 


f?7+0.32 
-0.08 


o.8ii°:«t 


0. 


00+0.18 
•'-'^-0.06 


91 


3197.5l^J^:^ 


0. 


1 Q+0.09 
•-'-°-0.06 





.0401°:°^ 


1 9Q+1.86 


0. 


91+0.16 
■''-^-0.08 


2 70-1-1-00 

Z. /U_2 7Q 


0. 


99+0.09 
•^^-0.06 


92 


4350.0lf,°:° 


2. 


qc:+0.90 
•■^^-0.68 


1, 


.1801°:!? 


2.42li:°^ 


1. 


cq+0.39 
■^''-0.45 


1 4Q+3-46 
1-^»_1.49 


1. 


19+0.54 
•1^-0.66 


93 


3270.0iy«:° 


0. 


1 1 +0.04 
•-^-^-0.03 


0, 


.0301°:°^ 


3-07li:^^ 


0. 


9f.+0.11 
■"^"-0.15 


^•^^-4.34 


0. 


90+0.06 
•^•^-0.10 



Table 4: Ha line widths and mass accretion rates. A positive EW value 
indicates emission. 



ID EW[Ha] CTTS? Hq [10%] CTTS? log Mac Double peak? 

(A) (km s-i) (Mo yr-i) 



7 


201.32 


yes 


351.0 ±5.0 


yes 


-9.49 ± 0.39 




13 


4.15 




98.0 ± 11.0 




-11.94 ±0.33 




16 


27.16 


yes 


306.0 ± 4.0 


yes 


-9.92 ± 0.37 


yes 


17 


4.66 




94.0 ±11.5 




-11.98 ±0.33 




18 


7.07 




116.0 ±7.5 




-11.76 ±0.32 




19 


2.23 




100.0 ± 11.0 




-11.92 ± 0.33 


yes 


20 


8.24 




137.0 ± 16.5 




-11.56 ±0.35 




22 


327.34 


yes 


312.0 ±4.5 


yes 


-9.86 ± 0.37 


yes 


23 


5.20 




183.0 ± 16.0 




-11.11 ±0.36 


yes 


24 


10.62 




115.0 ±6.0 




-11.77 ±0.32 




25 


7.01 




112.0 ± 7.0 




-11.80 ±0.32 




27 


36.56 




137.0 ± 5.5 




-11.56 ± 0.32 




47 


9.77 




101.0 ± S.O 




-ii.cS8 ± 0.32 
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YSOs in the Lupus Molecular Clouds 



Table 4: Ha line widths and mass accretion rates. A positive EW value 
indicates emission. 



ID 


EW[Ha] 

(A) 


CTTS? 


Ha [10%] 
(km s-i) 


CTTS? 


log Mac 

(Mo yr-i) 


Double peak? 


48 


3.37 




159.0 ± 7.0 




-11.35 ±0.33 




49 


44.74 


yes 


224.0 ± 19.5 


yes 


-10.72 ±0.39 


yes 


50 


147.24 


yes 


385.0 ±6.5 


yes 


-9.16 ±0.41 


yes 


51 


11.55 




125.0 ±6.0 




-11.68 ±0.32 




52 


130.76 


yes 


289.0 ±4.0 


yes 


-10.09 ±0.36 




54 


26.84 


yes 


501.0 ± 6.0 


yes 


-8.03 ± 0.47 


yes 


58 


10.66 




139.0 ±11.0 




-11.54 ±0.33 


yes 


59 


255.62 


yes 


264.0 ±8.0 


yes 


-10.33 ±0.36 




62 


22.29 


yes 


431.0 ± 10.5 


yes 


-8.71 ± 0.44 


yes 


65 


17.68 




463.0 ± 27.5 


yes 


-8.40 ± 0.52 


yes 


66 


23.59 




181.0 ± 6.0 




-11.13 ±0.33 




67 


9.43 




149.0 ± 6.0 




-11.44 ±0.32 




70 


11.93 




123.0 ±6.5 




-11.70 ±0.32 




71 


0.68 




167.0 ± 24.0 




-11.27 ±0.40 




76 


97.79 


yes 


348.0 ± 4.5 


yes 


-9.51 ± 0.39 




77 


23.74 


yes 


504.0 ± 9.0 


yes 


-8.00 ± 0.47 




78 


333.58 


yes 


283.0 ± 7.0 


yes 


-10.14 ±0.37 




79 


6.37 




171.0 ±9.5 




-11.23 ±0.34 




80 


26.84 


yes 


182.0 ±4.5 




-11.12 ±0.33 




81 


10.54 




236.0 ± 6.0 


yes 


-10.60 ±0.35 




82 


68.76 


yes 


389.0 ± 6.0 


yes 


-9.12 ±0.41 


yes 


83 


71.81 


yes 


440.0 ± 5.5 


yes 


-8.62 ±0.43 




84 


267.22 


yes 


450.0 ± 6.5 


yes 


-8.52 ± 0.44 


yes 


85 


139.22 


yes 


261.0 ±3.5 


yes 


-10.36 ±0.35 




86 


40.99 


yes 


486.0 ± 6.0 


yes 


-8.18 ±0.46 


yes 


87 


58.32 


yes 


229.0 ±5.0 


yes 


-10.67 ±0.34 




88 


121.73 


yes 


498.0 ± 4.0 


yes 


-8.06 ± 0.46 


yes 


89 


7.49 




156.0 ±6.0 




-11.38 ±0.32 




90 


11.40 


yes 


344.0 ± 6.5 


yes 


-9.55 ± 0.39 




91 


38.67 


yes 


395.0 ± 6.5 


yes 


-9.06 ± 0.41 




92 


26.42 


yes 


432.0 ± 4.0 


yes 


-8.70 ± 0.43 


yes 


93 


116.96 


yes 


379.0 ± 3.0 


yes 


-9.21 ± 0.40 





